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Abstract 

If the energy momentum tensor contains bulk viscous stresses violating the dominant energy 
condition (DOC) the energy spectra of the relic gravitons (produced at the time of the 
DOC's violation) increase in frequency in a calculable way. In a general relativistic context 
we give examples where the DOC is only violated for a limited amount of time after which 
the ordinary (radiation dominated) evolution takes place. We connect our discussion to some 
recent remarks of Grishchuk concerning the detectability of the stochastic gravitational wave 
background by the forthcoming interferometric detectors. 
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Every transition of the background geometry leads, inevitably, to the production of 
stochastically distributed relic gravitons whose energy spectra represent a crucial probe of 
the very early stages of the evolution of our Universe [|l]. Prior to the nucleosynthesis epoch 
there are no direct test of the thermodynamical state of the Universe and the presence of 
an inflationary phase of expansion is usually justified by causality arguments applied to the 
Cosmic Microwave Background (CMB) photons whose emission regions could not have been 
in causal contact in the far past if a never ending radiation dominated phase would precede 
our present matter dominated stage of expansion M . 

An inflationary evolution, if regarded at an effective level, violates the (general relativis- 
tic) strong energy condition (SEC) namely p + 3p > where p and p are, respectively, the 
energy density and the pressure density of the perfect fluid sources driving the evolution 
of the background geometry ||. The need for such a violation can be immediately seen 
from the structure of Friedmann equations which imply that, if a > 0, the SEC needs to be 
violated (a(t) is the scale factor of the homogeneous and isotropic Universe and the over-dot 
denotes the differentiation with respect to the cosmic time t). 

Recently, Grishchuk 0, ^ made an interesting observation concerning the slopes of the 
energy spectra of relic gravitons emerging from the models of the early Universe. In short 
the argument goes as follows. Suppose that we are in the framework of general relativity and 
suppose, as it is usually done, that the effective sources driving the background geometry can 
be parameterized by a stress tensor with perfect fluid form. As we stressed at the beginning, 
there are little doubts that the Universe had to be dominated by radiation (at least) since 
nucleosynthesis. Prior to that epoch we do not know which kind of energy momentum tensor 
could approximate the background sources but we would like to deal with ever expanding 
Universes. Moreover, at least for some time, we would like an accelerated expansion in order 
to solve the so called kinematical problems of the standard cosmological model. The second 
of the two previous requirements necessarily leads to the violation of the SEC. 

A very naive model of the early Universe would then be given by two phases. An unknown 
phase where the perfect fluid sources have a generic equation of state p = 7p , followed, at 
some transition time ti, by a radiation dominated phase with p = p/3. The question we 
are very interested to ask is under which conditions such a toy model would produce energy 
spectra of relic gravitons increasing faster than the first power of the frequency. Such a 
question is of obvious experimental relevance since, in the near future, various interferometric 
detectors of gravitational waves will come in operation. Now, if the energy spectra of relic 
gravitons are flat (or slowly increasing) with frequency there are little hopes of detecting 
them. In fact, the COBE limit applied to a frequency uq ~ 1CT 18 ho Hz imposes that the 
relic gravitons energy density (in critical units) /iq^qw has to be smaller then (or of the order 
of) 6.9 x 10~ n (ho is the present indetermination in the value of the Hubble constant). Due 
to the transition from radiation to matter, the infra-red branch of the graviton spectrum 
declines (in frequency) as v~ 2 between z/ and u 2 ~ 1CU 16 Qq Hz. If we take into account 
that the typical frequency of operation of the interferometers is between 10 and 100 Hz, 
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we see that it is quite interesting to understand which kind of models could give energy 
densities (around 100 Hz) larger than the typical inflationary models whose prediction is 
^0 GW < lO" 14 . 

If we confine our attention to expanding Universes (i.e. a > 0) the slopes of the energy 
spectra are crucially related to the sign of p+p. This point can be easily seen by looking at the 
Friedmann equations whose solution gives a(t) = t a with a = 2/[3(7 + 1)], having assumed 
p = 7p. The spectra of relic gravitons produced because of transition of the geometry from 
a generic perfect fluid stage to a radiation stage can be easily computed by matching the 
solutions of Eq. flllD in the two temporal regions defined by the transition time t%. The 
requirement that we want expanding (and inflationary) Universes implies that the energy 
density of the created gravitons cannot increase with frequency if p + p > 0. 

It was noticed by Barrow |J that a violation of the DOC does not necessarily forbid the 
viability of a given cosmological model in the same way as the violation of the SEC does not 
forbid the viability of ordinary inflationary models M. A violation of the DOC implies that 
H > since, by combining the Friedmann equations we get that 

M 2 P H=~( P + p). (1) 

Notice, moreover, that the requirement of eternal expansion implies, together with H > 0, 
that a > 0. In the context of general relativity, cosmological solutions violating the DOC 
(but not violating the weak energy conditions) have been sometimes connected with the 
presence of bulk viscous stresses driving the evolution of the geometry. 

Bulk viscous stresses || can appear quite naturally in cosmology as the result of the pro- 
cesses of particles production [|H|] in curved space-time . Moreover, bulk viscosity represents 
an interesting extension of the standard cosmological scenarios since it introduces dissipation 
in the model without breaking isotropy but by only redefining the space-space components 
of the energy momentum tensor according to 

p' = p- 3£H, (2) 

where £ is the bulk viscosity coefficient. The Friedmann equations in the presence of bulk 
viscous stresses can be easily written as 

M 2 P H 2 = p, 

Ml(H + H 2 ) = - l -(p + ?>pi), 

p + 3H(p + p') = 0. (3) 

As an example of the possibility of violating the DOC in not unreasonable models let us 
assume the simplest possibility namely £ = const > and 7 = — 1 which represents a 
particular case of the solution described by Barrow. The scale factor associated with this 
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solution has a typical double exponential form || 

9 

a(t) = a exp [exp (-ft)] (4) 

and it describes an accelerated expansion with H > (we took Mp = 1). In the context of 
|J it was argued that this behavior may be connected with a violation of the second principle 
of thermodynamics generalized to curved backgrounds (of quasi-de Sitter and non de Sitter 
type 0). Other solutions of the type given by Eq. @) can be constructed by assuming, for 
instance, £ = ap m , m being an arbitrary power and a a proportionality constant. 

A distinctive feature of the class of solutions by Eq. (jD is that they describe a violation 
of the DOC which occurs all the time. Namely, we can show that H > holds over all the 
times where the solution is defined. 

Solutions can be constructed where the violation of the DOC is only realized for a finite 
amount of time. The phase leading to a violation of the DOC can then be smoothly con- 
nected with an ordinary evolution preserving both the DOC and the SEC and whose Hubble 
parameter declines in time. In order to have this behavior dynamically realized we have to 
require that the bulk viscosity coefficient will change sign at the moment connected with 
the transition from the phase violating the DOC to the phase where the DOC is restored. 
Let us then assume that £ = qH/H where q is a constant to be determined by consistency 
with the Friedmann equations. The reason for this parameterization is indeed simple. As 
we stressed, the amount of violation of the DOC is indeed proportional to H. Moreover, the 
compatibility with the Eqs. (|3]) requires that q = (2/9)M P . 

Under this ansatz, a class of solutions of Eqs. (13) is given by the scale factor 



a(t) = exp [5 arcsinh(t)] = ( t + ^t 2 + tjj , (5) 
whose associated Hubble parameter is given by 

H =^=„- (6) 



/t 2 + 1\ 

We can immediately see by taking the limit of a(t) for t — > ±oo we get that 

a ±OQ (t) - (±t) ±s (7) 

In this model the explicit evolution of the shear parameter and of the energy density can be 
obtained, according to Eqs. (||), as 

Notice that £(£) is positive for t — > — oo and it is negative for t —>■ +oo. Therefore, the DOC 
is violated in the cosmic time interval ] — oo,0] and it is restored for [0, +oo[. The energy 
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density and the curvature are regular and well defined for every t and, moreover, a > for 
] — oo, 0] whereas a < for [0, +oo[. 

If we want the phase where the DOC is violated to be smoothly connected to a radiation 
dominated evolution we are led to choose 5 = 1/2. With this choice we can see from the 
exact expressions of the scale factor that 



a(t) = [t+ ji* + t 2 j , (9) 

Notice that, in this case, the pressure density goes as p — > p/3 for t — > +oo. The examples 
we just discussed suggest a possible speculation. The initial phase where the DOC and the 
SEC are both violated could evolve towards a phase where the SEC and the DOC are both 
restored. 

If this is the case we can show that the obtained graviton spectra excited by this type of 
background might be of some phenomenological relevance. Suppose indeed that the Universe 
evolves through two different stages of expansion. In the first stage the DOC and the SEC 
will be both violated (i.e. a(t) ~ (— t)^ 1 / 2 ) In the second phase the DOC and the SEC are 
both restored (i.e. a(t) ~ t 1//2 ). If the transition from the regime p + p < to the regime 
p + p > occurs smoothly we can also expect the presence of a transition epoch where 
p + p = and the DOC is saturated but not violated. Then by ensuring the continuity of 
the scale factors (and of their first derivatives) between the three regimes we have that the 
background evolution can be written, in conformal time, as 

a (v) = [ ] 3 7 V<-Vh ^p + p<0 

Vi 

oi = K > ~Vi<V<V2, ^p + p = 

3771 

The evolution equation of the proper amplitude of each polarization of the tensor fluctu- 
ations hij of a conformally flat background metric of Friedmann- Robert son- Walker type can 
be written, in conformal time and in Fourier space, as 

ip" + [k 2 - = 0, (11) 

CL 

where for each polarization if) = ah (the prime denotes derivation with respect to conformal 
time). 

The amplification induced by the background of Eq. (|T0|) in the proper amplitude of the 
gravitational waves can be obtained from the solutions of Eq. ( p]) in the three temporal 
regions separated by 771 and r\2- 

vr ,2^ 1 „( 2 ) 
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Mkv) = ^^[e^A+WH^ikv) + A-.{k)e*"H$(kn)], - Vi < V < i», 
1> nj (kn) = -}=[B + (k)e- ikr > + B^{k)J k \ V>V2 (12) 



where Hj) 1 ' 2 ^ are the Hankel functions. By matching in 77 = —771 ipijijn and ^ /7 m in — r7 X 
and 772 we can determine the mixing coefficient A±(k), B±(k) whose frequency behavior will 
fix the spectral evolution of the energy density. 

This is a very straightforward calculation |Tl| and the final result can be expressed in 



terms of the energy density of the relic gravitons at the present time and in critical units 

Qgw(^)*o) = ^ 7 (io)s 2 , Co>2 < uj < Ci?i 

4 

2 u \ 3 



tt GW (u,t ) = fi 7 (t )s , Wdec < ^ < ^2 (13) 

\UJ2/ 

where s ~ Hi/Mp; fi 7 ~ 2.6 x 10 -5 /ig 2 is the present fraction of critical energy density 
stored in radiation and u = k/a = lisv. 

Now our argument is very simple. The only free parameter in this model is given by 
the duration of the intermediate phase saturating the DOC, or, in other words by the ratio 
it)\ju)2- The overall normalization of the spectrum can be fixed by assuming, at the 
maximal amplitude compatible with the nucleosynthesis bound f2Gw~10 -5 . Taking into 
account that uj\ ~ 10 11 Hz we must conclude that the energy density of the relic gravitons 
at the scale of the interferometric detectors (i.e. u)i ~ 0.1 kHz) is given by 

n GW ~ 1(T 5 (^) 4 / 3 . (14) 
^2 

Therefore, if the duration of the intermediate phase lies in the range \r\%jr\]\ ~ 10 7 -10 8 we 
can conclude that the energy density of relic gravitons can be as large as 10~ 7 -10~ 6 in the 
frequency range probed by the interferometers without conflicting with other bounds coming 
from lower frequencies. 

We showed that the existence of tilted spectra of relic gravitons can be connected, in 
the context of general relativity, with the violation of the DOC induced by bulk viscous 
stresses. In many respects what we presented in this paper are toy examples which are 
only suggestive of a possible dynamics. In the near future these speculations can have some 
phenomenological impact. Possible indications of the existence of tilted gravitons spectra at 
the frequency of the interferometric devices can put bounds on the possible violation of the 
DOC in the early Universe. If tilted gravitons spectra of the type discussed in the present 
paper will be detected, then the conclusion is twofold. Either the DOC was violated in the 
early Universe, or the theory of gravity underlying the correct model describing the early 



stages of the evolution of the Universe was not of Einstein-Hilbert type [12 
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